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ABSTRACT 

High-resolution observations provide evidence about the existence of small- 
amplitude transverse oscillations in solar filament fine structures. These oscil- 
lations are believed to represent fast magnetohydrodynamic (MHD) waves and 
the disturbances are seen to be damped in short timescales of the order of 1 
to 4 periods. In this Letter we propose that, due to the highly inhomogeneous 
nature of the filament plasma at the fine structure spatial scale, the phenomenon 
of resonant absorption is likely to operate in the temporal attenuation of fast 
MHD oscillations. By considering transverse inhomogeneity in a straight flux 
tube model we find that, for density inhomogeneities typical of filament threads, 
the decay times are of a few oscillatory periods only. 

Subject headings: MHD — Sun: filaments — Sun: oscillations — waves 



1. Introduction 

Quiescent solar filaments form along the inversion p olarity line or between the weak 
remnants of active regions. Early filament observations (jEngvoldl Il998l ) as well as recent 
high-re solution Ha ob servations, obtained with the Swedish Solar Telescope (SST) in La 
Palma (ILin et al.ll2005l ) and the Dutch Open Telescope (DOT), have revealed that their fine 
structure is composed by many horizontal and thin dark threads. The measured average 
width of resolved threads is about 0.3 arc sec (~ 210 km) while the length is between 5 
and 40 arc sec (~ 3500 - 28,000 km). They seem to be partially filled with cold plasma 



( ILin et al.ll2005l ). typically two orders of magnitude denser than t hat of the corona, and it 



is generally assumed that they out line their magnetic flux tubes (jEngvoldl Il998l : iLinl 12004 



Lin et al.l 120051 ; iMartin et al.ll2008l ). This idea is strongly supported by the fact that they 
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are inclined with respe ct to the filament long axis at a similar angle to what has b een found 
for the magnetic field (lLeroylll980l ; iBommier et al.lll994t iBommier fc Leroyl 119981 ) . 



Small amplitude oscillations have been observed in filaments and it is well established 
that these periodic changes are of local nature. The detected peak velocity ranges from 
the noise level (down to 0.1 km s _1 in some cases) to 2-3 km s _1 , a lthough larger values 



have also been repor ted. Two-dimensional observations of filaments by lYi &: Engvoldl (1l99ll ) 
and lYi et al.l (I199ll ) revealed that individual threads or groups of threads may oscillate 
independe ntly with their own periods, which range between 3 and 20 minutes. More recently, 
Linl (120041 ) reports that spatially coherent oscillations are found over slices, with an area of 
1.4 x 54 arc sec 2 , of a polar crown filament, and that among other, a si gnificant periodic ity at 
26 minutes, strongly damped after 4 periods, appears. Furthermore, iLin et al.l (120071 ) have 
shown evidence about traveling waves along a number of filament threads with an average 
phase velocity of 12 km s _1 , a wavelength of 4" (~ 2800 km), and oscillatory periods 
of the individual threads that vary from 3 to 9 minutes. The observed periodic signals 
are obtained from Doppler velocity measurements and can therefore be associated to the 
transverse displacement of the fine structures. 

The observed small amplit ude oscillations have be en interpreted in terms of magne- 
tohydrodynam ic (MHD) waves (jOliver &: Ballesterl 120021) and theo r etical models have been 
developed (see Ballester 2005 . 20061 . for recent reviews) . Diaz et al. ( 2002 ) modeled a promi- 
nence thread as a straight cylindrical flux tube with a cool region representing the filament 
material, confined by two symmetric hot regions. They found that the fundamental fast 
mode is always confined in the dense part of the flux tube, hence, for an oscillating cylindri- 
cal filament thread, it should be difficult to induce oscillations in adjacent threads, unless 
they are very close. 

The time damping of prominence oscillations has been unambiguously determined in 
some observations. Reliable val ues for the damping time ha ve be en derived, from diffe rent 
Doppler velocity tim e ser i es, by iMolowny-Horas et al.l (Il999l ) and iTerradas et al.l (120021 ). in 
prominences, and by iLinl (120041 ). in filaments. The values thus obtained are usually between 
1 and 4 times the corresponding period, and large regions of prominences/filaments display 
similar damping times. The damping of perturbations is probably a common feature of fil- 
ament oscillations, hence theoretical mechanisms must be explored and their damping time 
scales should be compared with those obtained from observations. Linear non-adiabatic 
MHD waves have been proposed as a potential mechanism to explain the observed attenu - 
ation time scales JCarbonell et al.ll2004j ; ITerradas etaDboOlL I2OO5I : ISoler et al.lbo07l . l2Q08h . 
Using thermal mechanisms only slow waves can be damped in an efficient manner while 
fast waves remain almost undamped. Ion-neutral collisions provide a possible mechanism to 



- 3- 



damp fast waves (as well as Alfven waves) that is able to reprodu ce observed dampin g times 
for given parameter values, in particular for a quasi-neutral gas (IForteza et al.l 120071 ). 



Apart from the mentioned non-ideal damping mechanisms, there is another possibility to 
attenuate fast waves in thin filament threads. The phenomen on of resonant wave damping is 
well d ocumented for fast kink waves in coronal loops (see e.g. iGoossens et al.ll2006l ; iGoossens 
20081 . for recent reviews) and provides a plausible explanation t o quickly damped transv erse 
loop oscillations observed by TRACE (INakariakov et al.l 11999c lAschwanden et al.lll999l ). In 
this Letter, we address the resonant damping mechanism in the context of filament thread 
oscillations and assess its relevance in explaining the observed attenuation time scales. 



2. One-Dimensional Non-Uniform Filament Thread Model 



Given the relatively simple structure of filament threads, when compared to the full 
prominence/filament structure, the magnetic and plasma configuration of an individual and 
isolated thread can be theoretically approximated using a rather simplified model. We 
consider a gravity-free, straight, cylindrically symmetric flux tube of mean radius a (see 
Fig- 1])- In a system of cylindrical coordinates (r, <p, z) with the z-axis coinciding with the 
axis of the tube, the magnetic field is pointing in the ^-direction, B = Be z . We neglect 
gas pressure, which allows us to concentrate on the oscillatory properties of fast and Alfven 
MHD waves and their mutual interaction. In our straight field configuration this zero-/? 
approximation implies that the field strength is uniform and that the density profile can 
be chosen arbitrarily. The inhomogeneous filament thread is then modeled as a density 
enhancement with a one- dimensional non-uniform distribution of density, p(r), across the 
structure. The internal filament plasma, with uniform density, pf, occupies the full length 
of the tube and is connected to the coronal medium, with uniform density, p c , by means of 
a non- uniform transitional layer of thickness I. The ratio l/a provides us with a measure of 
the transverse inhomogeneity length-scale, that can vary in between l/a = (homogeneous 
tube) and l/a = 2 (fully non- uniform tube). 



3. Damping of Linear Fast Kink Waves 



The magnetic flux tube model adopted here is a waveguide for a number of MHD 
oscillatory solutions. In the zero-/? approximation slow waves are absent. The properties of 
the remaining small amplitude fast and Alfven waves can readily be described by considering 
the linear MHD wave equations for adiabatic changes of state for perturbations of the form 
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f (r)exp(i(u)t + rmp — k z z)). Here m and k z are the azimuthal and longitudinal wave-numbers 
and 10 the oscillatory frequency. We further concentrate on perturbations with m — 1, 
which represent fast kink waves that produce the transverse displacement of the tube as 
they propagate along the density enhancement. 



3.1. Analytical theory 



In the long wavelength or thin tube approximation ( k z a <C 1) the frequency of the 
m = 1 fast kink wave can be written down analytically (see lEdwin fc Robertslll983l ) as 



Pf + Pc 



(i; 



with Va/,c = B/ y/HPf, c the filament (f) and coronal (c) Alfven velocities. By defining c = 
Pf/pc, for the density contrast, the period of kink oscillations with a wavelength A = 2n/k z 
can be written as 



y/2 A 



2 V, 



Af 



1/2 



(2) 



The factor containing the density contrast varies between \/2 and 1, when c is allowed to vary 
between a value slightly larger that 1 (extremely tenuous thread) and c — > oo. For typical 
filament thread densities, which are two orders of magnitude larger than coronal densities, 
this factor is near unity. Equation (JSJ) then predicts Alfven velocities in the thread as low 
as a few km s~ x or as large as ~ 200 km s _1 , for combinations of periods and wavelengths 
between 3 and 20 minutes, and 3000 and 20, 000 km, respectively. 

For fast kink waves to be damped by resonant absorption, transverse inhomogeneity 
in the Alfven velocity has to be considered. In our uniform field model this is obtained by 
considering I ^ 0. Then, the m = 1 solution is resonantly coupled to local Alfven waves. 
The coupling produces the temporal attenuation of fast transverse motions which are con- 
verted into localized azimuthal Alfvenic oscillations. Asymptotic analytical expressions for 
the damping time, r^, can be obtained under the assumption that the transverse inhomo- 
geneity length-scale is small (I /a 1). This is the so-called thin boundary approximation. 



written as (see e.g. 
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Td 
P 



a c + 1 

T c- r 



(3) 



Here F is a numerical factor that depends on th e particular variation of the density in the non - 



uniform layer. For a linear variation F = 4/-7T 2 (IHollweg &: Yang 19881 ; iGoossens et al.lll992l ); 
for a sinusoidal variation F = 2/n ( Ruderman fc Roberts! 2002 ). For example, considering 
c = 200 a typical density contrast and I /a = 0.1 equation (J3j predicts a damping time of ~ 
6 times the oscillatory period. 

Figure [2] shows analytical estimates computed using equation ([3]) (solid lines). The 
damping is affected by the density contrast in the low contrast regime and r^/P rapidly 
decreases for increasing thread density (Fig. [2^). Interestingly, it stops being dependent 
on this parameter in the large contrast regime, typical of filament threads. The damping 
time over period is independent of the wavelength of perturbations (Fig. [2b), but rapidly 
decreases with increasing inhomogeneity length-scale (Fig. [2b). These results suggest that 
resonant absorption is a very efficient mechanism for the attenuation of fast waves in filament 
threads, especially because large thread densities and transverse plasma inhomogeneities can 
be combined together. 



3.2. Numerical results 



The two approximations used to derive equations (|2J) and ([3]) may impose limitations 
to the applicability of the obtained results to filament thread oscillations. To start with, 
it is not clear how accurate the long wavelength approximation can be, especially for short 
wavelengths. The period in equation ([3]) corresponds to the long wavelength limit d2J), which 
is independent of the radius of the structure, and does not include effects due to radial 
density inhomogeneity. Also, equation (j3J) is only valid as long as the damping time is 
sufficiently larger than the period, an assumption made in order to derive the expression, 
and clearly cont radicted by the obtaine d resul ts. The accuracy of the damping formula (J3j 
was assessed by IVan Doorsselaere et al.l (120041 ) in the relatively low density contrast regime 
corresponding to damped coronal loop oscillations. Some degree of inaccuracy in the large 
contrast regime characteristic of filament threads is also to be expected. We have therefore 
computed numerical approximations to the solutions by solving the full set of linear, resistive, 
small amplit ude MHD wave equa tions for the m = 1 trans verse kink o scillations (see e.g. 



eqs. [l]-[5] in lTerradas et all 120061 ). using the PDE2D code (ISewelll 120051 ) . 



Figure [2] shows the obtained numerical results. Analytical and numerical solutions 
display the same qualitative behavior with density contrast and transverse inhomogeneity 
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length-scale (Figs. [2^, c). Now the damping time over period slightly depends on the wave- 
length of perturbations (Fig. [2b). Equation (j3J) underestimates/overestimates this magnitude 
for short/long wavelengths. The differences are small, of the order of 3% for c = 10, and 
do not vary much with density contrast for long wavelengths (A = 200a), but increase until 
6% for short ones (A = 30a). The long wavelength approximation is responsible for the 
discrepancies obtained for thin non- uniform layers (Fig. [2b) . Figure [2U shows how accurate 
equation ([3]) is for different combinations of wavelength, density contrast, and inhomogeneity 
length-scale. For thin layers (I /a = 0.1) the inaccuracy of the long wavelength approxima- 
tion produces differences up to ~ 10% for the combination of short wavelength with high 
contr ast thread. For thick layers, d ifferences of the order of 20% are obtained (in agreement 
with IVan Doorsselaere et al.ll2004l ). Here, the combination of large wavelength with high 
contrast thread produces the largest discrepancy. Numerical results allow the computation 
of more accurate values, but do no change our previous conclusions regarding the efficiency 
and properties of resonant damping of transverse oscillations in filament threads. 



4. Discussion 

In this Letter we have shown that due to the highly inhomogeneous nature of filaments at 
their transverse scales the process of resonant absorption is an efficient damping mechanism 
for fast MHD oscillations propagating in these structures. The relevance of the mechanism 
has been assessed in a flux tube model, with the inclusion of transverse inhomogeneity in the 
Alfven velocity. Also, the accuracy of analytical estimates, in terms of wavelength, density 
contrast, and transverse inhomogeneity, has been quantified. For the typical large filament 
to coronal density contrast the mechanism produces rapid damping in timescales of the order 
of a few oscillatory periods only. The obtained damping rates are only slightly dependent 
on the wavelength of perturbations. An important result is that the damping rate becomes 
independent of density contrast for large values of this parameter. This has two seismological 
consequences. First, the observational determination of density contrast is less critical than 
in the low contrast regime. Second, according to sei smic inversion resul t s that combine 



theor etical and observed periods and damping times (lArregui et al.l 120071 ; iGoossens et al. 



20081 ). high density thread models would be compatible with relatively short transverse 
inhomogeneity length-scales. Analytical estimates of l/a ~ 0.15 can even be calculated 
using equation for a given observed r d /P = 4, taking the limit c — > oo. 

We believe our conclusion on the relevance of resonant damping in filament thread 
oscillations is robust in front of the main simplification adopted for the present study. Ob- 
servations show that threads are only partially filled with cool and dense plasma and our 
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one-dimensional model misses this property. Although the resonant damping mechanism 
relies on the Alfven velocity inhomogeneity in the direction transverse to the magnetic field, 
an investigation of the damping properties in such a two-dimensional configuration is a nec- 
essary issue. Finally, the presence of flows is commonly observed along filament threads and 
the interplay between resonant damping and flows must also be explored in this context. 

The authors acknowledge the funding provided under projects AYA2006-07637 (Span- 
ish Ministerio de Education y Ciencia) and PCTIB2005GC3-03 (Conselleria d'Economia, 
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B 

Fig. 1. — Sketch of the model used to represent a radially non- uniform filament fine structure 
of mean radius a, with a transverse inhomogeneity length-scale I (see detail in the cross 
section). 
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Fig. 2. — (a)-(c): Damping time over period for fast kink waves in filament threads with 
a = 100 km. In all plots solid lines correspond to analytical solutions given by equation (J3]), 
with F = 2/tt. (a): As a function of density contrast, with I /a = 0.2 and for two wavelengths. 
(b): As a function of wavelength, with l/a = 0.2, and for two density contrasts, (c): As a 
function of transverse inhomogeneity length-scale, for two combinations of wavelength and 
density contrast, (d) Percentage difference, A, with respect to analytical formula ([3]) for 
different combinations of wavelength, A = 30a (dashed lines); A = 200a (dash-dotted lines), 
and density contrast. 



